Modulation of breathers in cigar-shaped Bose-Einstein condensates 



W. B. Cardoso, 1 A. T. Avelar, 1 and D. Bazeia 2 

1 Institute) de Fisica, Universidade Federal de Golds, 74.001-970, Goidnia, Golds, Brazil. 
2 Departamento de Fisica, Universidade Federal da Paraiba, 58.059-900, Joao-Pessoa, Paraiba, Brazil 

We present new solutions to the nonautonomous nonlinear Schrodinger equation that may be realized through 
convenient manipulation of Bose-Einstein condensates. The procedure is based on the modulation of breathers 
through an analytical study of the one-dimensional Gross-Pitaevskii equation, which is known to offer a good 
theoretical model to describe quasi-one-dimensional cigar- shaped condensates. Using a specific Ansatz, we 
transform the nonautonomous nonlinear equation into an autonomous one, which engenders composed states 
corresponding to solutions localized in space, with an oscillating behavior in time. Numerical simulations 
confirm stability of the breathers against random perturbation on the input profile of the solutions. 
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I. INTRODUCTION 

Breathers are solutions of nonlinear equations with a pro- 
file which is usually localized in space and periodic in time, 
even for equations with constant coefficients. Such solutions 
usually come from a compound of two or more solitons lo- 
cated at the same spatial position. They have been origi- 
nally introduced in the sine-Gordon equation [1], but they 
can also be found in other scenarios, controlled by the mod- 
ified Korteweg-de Vries equation [2], the Davey-Stewartson 
0] and the nonlinear Schrodinger equation Q]. In various 
physical systems, which are well described by nonlinear equa- 
tions, breathers directly affect their electronic, magnetic, opti- 
cal, vibrational and transport properties, such as in Josephson 
superconducting junctions |5, 6], charge density wave systems 
01, 4-methyl-pyridine crystals metallic nanoparticles 
conj ugated polymers 1 10], micromechanical oscillator arrays 
ifTTn . antiferromagnetic Heisenberg chains (13, [3, and semi- 
conductor quantum wells fllll . 

After the first experimental manipulation of Bose-Einstein 
condensates (BECs) fl5l - [l7ll . studies on bright and dark soli- 
tons fl8l - [2lll have triggered a lot of new investigations, with 
a diversity of scenarios being proposed and tested 1141 12214261 . 
In particular, the presence of experimental techniques for ma- 
nipulating the strength of the effective interaction between 
trapped atoms 1 27] leads us to believe that in BECs we have an 
excellent opportunity to investigate breathers of atomic mat- 
ter waves taking advantage of Feshbach-resonance manage- 
ment [28-36]. Thus, in this work the main aim is to show that 
breathers can be modulated in BECs, presenting a diversity of 
patterns. We offer an analytical study of such breather solu- 
tions of the Gross-Pitaevskii equation (GPE) with both the po- 
tential and the parameters that control the strength of the non- 
linearities being time-dependent. This makes the GPE nonau- 
tonomous and very hard to solve in general. However, we 
employ an Ansatz which transforms the nonautonomous GPE 
into an autonomous nonlinear Schrodinger equation (NLSE) 
which is easier to be studied. 

Since the one-dimensional GPE is a very good theoretical 
model to describe quasi-one-dimensional cigar-shaped BECs 
, we then investigate the presence of the breather be- 
havior in BECs under the assumption that a s <C a±_ <C A <C 



a 1 1 , for a s being the scattering length, A the spatial scale of the 
wave packet, and a± and a\\ being the characteristics trans- 
verse and longitudinal trap lengths, respectively. Also, in this 
work we use potentials and nonlinearities which are typical 
of BECs, and we believe that the results obtained below may 
stimulate new experiments in condensates. 

In the recent literature, attention has also been devoted to 
the case of discrete breathers, in arrays of BECs ll37l - [39ll . 
However, when one looks for continuous dynamics, with the 
GPE having variable coefficients depending on the space and 
time coordinates, analytical solution is not a trivial issue. 
Here we can point a numerical study ll34ll which emphasizes 
the stability properties of the breather excitations in a three- 
dimensional BEC with Feshbach-resonance management of 
the scattering length and confined by a one-dimensional opti- 
cal lattice. We believe that such excitations are made localized 
by the addition of an external potential, usually created by an 
optical lattice used to trap them in the lattice. 

In the present paper, however, we use a genuine breather 
solution, i.e., with oscillating behavior in time even when the 
nonlinear equation presents only constant coefficients (i.e., 
without modulation). The profile of these breathers are shown 
in Fig. [T] We also show that they can be modulated in space 
and time, controlled by the trapping potential and the strength 
of the cubic nonlinearity, as we show below for a diversity of 
situations. The procedure is similar to the former investiga- 
tions IN{)l445ll , but here we allow for the autonomous NLSE to 
be partial differential equation, and this introduces an interest- 
ing difference which we explore below. 

We start the investigation in the next section, where we 
present the general formalism, following as close as possi- 
ble the two recent works HHHl]. In Sec. Ill we illustrate the 
formalism with distinct examples, identified by the trapping 
potential and cubic nonlinearity, with a diversity of solutions. 
We also show in Sec. IV, through numerical simulation, that 
the solutions found in Sec. Ill are all stable under random per- 
turbation of 5% in their input profile. There we also com- 
ment of the case of nonpolynomial nonlinearity, which can be 
used to compensate for the one-dimensional approximation of 
real condensates. There we show that the one-dimensional ap- 
proximation used in this work is reliable. We finish the work 
in Sec. V, where we introduce the ending comments. 
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II. FORMALISM 

We start with the one-dimensional GPE with the general 
form 

#i = -\l>xx + V(X, t)l/> + P(|Vf)V>, CD 

where t) is the potential and P(\i(j\ 2 ) is given by 

TV 



wi 2 ) = E^+i(*)M 2n - 



(2) 



Thus, using ® in d6ct one gets p = yfa. Next, we rewrite the 
potential in the form 



V(x,t) = h{t)x 2 + f 2 (t)x + h{t), (10) 



with 



2a t 6 t 



(Ha) 



(lib) 



The coefficients g2n+i(t) of the polynomial function above 
describe the strength of the nonlinearities present in the sys- 
tem. Here we are using standard notation, in which x, t and 
all other quantities are dimensionless. 

Our goal is to find breather solutions of the Eq.©. Using 
the Ansatz 



4(x,t) = p(t)e i "( !E ' t )#(C(a:,t) ) r(t)) ) 



(3) 



we can change the nonautonomous GPE for an autonomous 
NLSE that may engender the breather solutions; it is given by 



i$ T = -I$ fc + Q(|*| 2 )*, 



where Q(|<£| 2 ) has the form 



N 



Q(|$| 2 ) = ^G 2n+1 |$ 



2n 



(4) 



(5) 



The coefficients G<m+\ are now constant parameters. Using 
(0 in ((TJ leads to ©, for p, rj, ( and r obeying the following 
equations 



T t = c 

Ct 



(6a) 
(6b) 

(6c) 



Here the potential and nonlinearities assume the form 

V(x,t) = -r h -r ] l (7) 

and 

fl 2i+ i(t) = ^^ ) i= 1,2,3,... (8) 

We introduce a(t) = ^Jr~ t to get from ([6ab the relation £ = 
a(t)x + 6(t). This is important because through f we can 
now determine the width of the localized solution in the form 
1 / a(t) and its center-of-mass position as —b(t) /a(t); see, e.g., 
Ref. [43]. In the same way, through the use of (l6bt one can 
write 



T](x, t) = —^X 2 — —X + c(t). 



(9) 



2a 2 ' 



(He) 



In this case, g2i+i{t) = G2i+ia 2 ~ % , for i = 1,2,3, and 
the nonlinearities assume a simple form. However, we note 
that in the above general procedure one needs that g$ = G5, 
suggesting that the strength of the quintic nonlinearity of the 
GPE (Q]) has to be constant. 

Similar ideas have been used recently in Refs. |4C 
There one changes the nonautonomous NLSE to an au- 
tonomous NLSE depending on a single coordinate, in the form 
of a static -like ordinary differential equation supporting soli- 
ton solutions. Such soliton solutions can be modulated via 
spatial and temporal dependences associated with the poten- 
tial and nonlinearity present in the nonautonomous NLSE. In 
this way, the soliton can be modulated in a time-oscillating 
pattern, as a breather solution; see, e.g., Ref S for a re- 
cent investigation. In our procedure, however, we maintain 
the nonlinear equation in the form given by Eq. ©, which is a 
partial differential equation which possesses breather solution 
even in the absence of the spatial and temporal modulations 
of the potential and cubic nonlinearity. This introduces an in- 
teresting difference because now we are directly controlling 
the breather solutions. We explore this possibility in the next 
section, investigating several distinct models. 



III. SPECIFIC MODELS 

Let us now turn attention to the case of cubic nonlinear- 
ity, with the strength gs(t) = g(t) and no other nonlinearity, 
which occurs for BECs with strong two-body scattering such 
as Rb or Na by use of Feshbach-resonance. In this case, the in- 
tegrability of the NLSE © for cubic nonlinearity allows that 
one obtains breather solutions with distinct profile. A spe- 
cific form of the two- soliton breather solution is obtained for 
G3 = G = — 1, which corresponds to the explicit solution 
El 



*(C,t) 



r/2 



4(cosh(3C) + 3e 4 * r cosh(C))e 
(cosh(4C) + 4 cosh(2C) + 3 cos(4r)) ' 



(12) 



It presents the simple form 3>(C,0) = 2sech(£) at r = 0. 
This is the simplest case, which we use to investigate distinct 
models below. 
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FIG. 1: (Color online) Breather solutions of the GPE ([T]), in the case 
of a vanishing potential and constant cubic strength. In (a) we display 
the case a = 1 and 6 = and in (b) a = 1 and b = t, making the 
breather to move in space. 



FIG. 3: (Color online) Plots of the (a) breather solution and (b) its 
profile, for the potential depicted in Fig. [2 
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FIG. 2: (Color online) Plot of the potential shown in JT3l ), using the 
parameters a — 1/(1 + cos 2 (t)) and 6 = 0. 



A. Vanishing potential 

Firstly, we study the free evolution with vanishing poten- 
tial and constant cubic nonlinearity strength. This case gives 
the standard breather solution, with the NLSE with zero po- 
tential and constant coefficients. Thus, we consider a = 1, 
which corresponds to g = —1 and j\ = 0, and we take 
ct = —b\ /2a 2 , leading to fs = 0. We can make f'2 = 
for b = and for b = t. These breathers are standard solu- 
tions and they are plotted in Figs. [T^ and[TJ), respectively. We 
note from (Q]) that for b = 0, the breather is localized in space 
and oscillates in time; however, for b = t the corresponding 
center of mass moves in space with oscillating time behavior. 
This behavior is related to Eq. ©, and to the nonzero phase in 
Eq. ©, together with the fact that here p is a constant and the 
form of the solution in Eq. (fT2l) also contributes with a non- 
trivial phase: the overall effect is that b = t makes the conden- 
sate to move away from its starting position, as it is explicitly 
shown in Fig.Q}). The oscillation period and frequency of the 
breather is of approximately 1.57 and 0.64, respectively. We 
obtain the limits for minimum and maximum amplitude equal 
to 4 and 16, respectively, and the time width at half height is 
approximately 0.39. Recall that we are using dimensionless 
coordinates. 



B. Flying bird potential 

Next, we deal with the specific case of a potential with 
space and time dependence, and a cubic nonlinearity with 
time-dependent strength. Here we choose the nonlinearity in 
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FIG. 4: (Color online) Plots of the (a) breather solution and (b) its 
profile. Here we are using a and b as in Fig. [2] 



the form # = —1/(1 + cos 2 (t)). To get a simple potential we 
take 6 = and c = 0, which make f'2 = and fs = 0. In this 
case we get to the potential 



2cos(2£) 2 

V — — - — x 

3 + cos(2t) 



(13) 



Figs. [2] and [3] show the potential and the breather solution, 
|^| 2 , respectively. This potential changes from attractive to 
expulsive behavior controlled periodically 1147144911 . We call 
this the flying bird potential, since it simulates the wings of a 
flying bird. Here we note that, differently from the first case, 
the periodic modulation does affect the periodic behavior of 
the breather solution. This interesting behavior appears from 
Figs.Q]and[3l there one sees that the breather exhibits a reduc- 
tion in its oscillation frequency, which is modulated through 
the trapping potential and the strength of the cubic nonlinear- 
ity. 

We see that the oscillation period and frequency of the 
breather is of approximately 3.14 and 0.32, respectively. Note 
that the breather frequency was reduced by half when com- 
pared with the former case. Also, we obtain the limits for 
minimum and maximum amplitude equal to 2 and 16, respec- 
tively, and the time width at half height is approximately 0.54. 

The procedure can induce another important behavior to 
the solution. We can make the breather to move with the 
choice of b such that b = - tan(t)/16(tan 2 (t) + 2) + 
3\/2 arctan(tan(t)/v^2)/32, to get fi = 0; also, the choice 
c t = —b 2 /2a 2 makes fs = 0. In this way, using the same non- 
linearity above the potential is also the same shown in (TT3l) . 
But now the breather present a nontrivial phase in ©, and it 
presents an interesting behavior: it seems to split in two parts 
in the time coordinate. This is different from the spatial split- 
ting which was studied before in ll35l l36h . Our solution is 
depicted in Fig.|4l and the breather seems to have two distinct 
time scales. 
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FIG. 5: (Color online) Plots of the (a) potential shown in dT4t and 
(b) the corresponding moving breather solution for a constant and 
b = sin(t). 
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FIG. 6: (Color online) Plot of the potential shown in i\5\ . 



C. Seesaw potential 

Another example is to verify the behavior of the moving 
breather solution in presence of a potential with linear spatial 
dependence. To this end we choose j\ = 0, through the use of 
constant nonlinearity a = 1, with fs = for c t = —b 2 /2a 2 . 
Thus, b = — sin(t) corresponds to the potential 



V = sm(t) x. 



(14) 



This potential has a seesaw behavior, as it is displayed in 
Fig. [5^. In Fig. 03 we show the zig-zag pattern of the cor- 
responding breather solution. In this case, the period of os- 
cillation and frequency of the breather solution are 1.57 and 
0.64, respectively; note that these values are the same of the 
first case. The amplitude now oscillates between 4 and 16 and 
the time width at half height is approximately 0.41. 



D. Another potential 

Other potentials can be studied for different values of a and 
b. For instance, using the cubic nonlinearity strength g = 
— 1/(1 + cos 2 (t)), with b = t, and making fs = we obtain 
the following potential 



V = 



2 cos(2£) 
3 + cos(2t) 



[1 + 5cos 2 (t)]sin(t)x. (15) 



It is displayed in Fig. [6] In this case, the breather solution 
presents an interesting moving pattern, as shown in Fig. [7^, 
for b = sin(t). Also, in Fig. [7J) we display the behavior of 
the moving breather considering b = s'm(y/2t). With these 
choices, the breather can present periodic or quasiperiodic 
pattern, respectively, as it appears explicitly in the Figs. [7^ 
and|7]D. In the case of periodic pattern, we obtain the oscilla- 
tion period of 3.14 and frequency of 0.32, with the time width 
at half height given by 0.56. 



IV. NUMERICAL SIMULATIONS 

Starting with the GPE (Q} in the case of cubic nonlinearity, 
we employ a numerical study focusing on the propagation of 
the modulated breather solution in the form given by Eq. ([3]). 
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FIG. 7: (Color online) Plots of the profile of the moving breather 
solution, showing the case of (a) periodic [b = sin(t)] and (b) 
quasiperiodic [b = sin(v / 2t)] patterns. 



The numerical method is based on the split- step finite differ- 
ence algorithm with a Crank-Nicolson method to solve the 
dispersive term in (Q]). It is known that the Crank-Nicolson 
method is unconditionally stable II 5 Oil. In this way, using a 
short length of time one can obtain the true pulse propaga- 
tion. Here, we used the time and space step-size of 0.0001 
and 0.01, respectively. Thus, with a random perturbation of 
5% in the input profile the solution can expose its stable or 
unstable behavior. All the solutions were tested and appeared 
to be stable under the suggested random perturbations. We 
analyzed the solutions for a very long time, until the rescaled 
time t = 10000. The results strongly suggest experimental 
verification of the above breather excitations in condensates. 

Another important question to be considered is related to 
the validity of the one-dimensional approximation to BECs, 
which we are using in this work. Evidently, we are assuming 
the regime a s <C a±_ <C A <C a\\, where the Eq. (Q]) is sup- 
posed to be valid. However, a more general possibility appears 
when one thinks of considering the quasi-one-dimensional dy- 
namics, which can be treated through a nonpolynomial non- 
linear equation, as shown, for instance, in the recent works 
ll5ll - l55ll . In this case, the Ansatz used in Eq. © does not 
transform the nonautonomous equation into an autonomous 
one. 

To get further information on this issue, let us compare the 
results obtained in the present study with a numerical inves- 
tigation of the quasi-one-dimensional equation in presence of 
nonpolynomial nonlinearity given by Dll455ll 



l + (3/2)g(M)M' 



(16) 



Note that for a weak-coupling regime, i.e., g(x, t)\ifj\ 2 <C 1, 
one can easily obtain the effective one-dimensional Eq. (Q]); 
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FIG. 8: (Color online) Plots of max \ip\ 2 for comparison between 
the cubic (orange, dashed line) and nonpolynomial (blue, solid line) 
nonlinearities, considering g = —0.1 and V = 0. 

see e.g., ll54ll for more details on this issue. In our model, this 
weak-coupling regime is contemplated when one decreases 
a(t), i.e., when one increases the width of the localized so- 
lution. As a simple example, let us compare the case A of a 
vanishing potential, studied above, with that obtained numer- 
ically through Eq. dT6b . We use the solution (Tf2l) at t = 
as initial input in dT6b with a = 0.1 and 6 = 0, which cor- 
respond to g — —0.1 and V = 0. The results are shown in 
Fig. [U Note that the solid (blue) curve presents a slight de- 
lay, compared to the dashed (orange) curve which represents 
case A, of vanishing potential. We believe that this delay is 
due to the max (gl^l 2 ) — 0.16, which should be interpreted 
as <C 1. The results depicted in Fig. [8] in solid (blue) line and 



in dashed (orange) line reinforces the present study, showing 
that the one-dimensional approximation used in this work is 
reliable. 



V. ENDING COMMENTS 

In summary, we have studied the modulation of breather ex- 
citations in a nonautonomous GPE. We have used an Ansatz to 
transform the GEP into an autonomous NLSE which supports 
breathers. The procedure allowed us to investigate a diversity 
of solutions, each one with its distinctive behavior. The stan- 
dard breather behavior appears in the simple case, with van- 
ishing potential and constant nonlinearity strength. But we 
could find other features, such as the moving, and the periodic 
or quasiperiodic behavior. 

We have studied stability numerically, with the results 
showing that all the breather excitations which we have found 
are stable against random perturbation on the input profile of 
the solutions. Also, we have investigated a simple case of 
nonpolynomial potential, and the results show that the one- 
dimensional approximation which we have used in the work is 
reliable. The present study opens several new possibilities, in 
particular on the experimental search for breather excitations 
in cigar- shaped condensates, and on extensions of the proce- 
dure to condensates with more general nonlinearities and in 
higher dimensions. 
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